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ABSTRACT 

There has been a persistent conundrum in attempts to model the nucleosynthesis of heavy elements 
by rapid neutron capture (the r-process). Although the location of the abundance peaks near nuclear 
mass numbers 130 and 195 identify an environment of rapid neutron capture near closed nuclear 
shells, the abundances of elements just above and below those peaks are often underproduced by 
more than an order of magnitude in model calculations. At the same time there is a debate in the 
literature as to what degree the r-process elements are produced in supernovae or the mergers of 
binary neutron stars. In this paper we propose a novel solution to both problems. We demonstrate 
that the underproduction of nuclides above and below the r-process peaks in main or weak r-process 
models (like magnetohydrodynamic jets or neutrino-driven winds in core-collapse supernovae) can 
be supplemented via fission fragment distributions from the recycling of material in a neutron-rich 
environment such as that encountered in neutron star mergers. In this paradigm, the abundance 
peaks themselves are well reproduced by a moderately neutron rich, main r-process environment 
such as that encountered in the magnetohydrodynamical jets in supernovae supplemented with a 
high-entropy, weakly neutron rich environment such as that encountered in the neutrino-driven-wind 
model to produce the lighter r-process isotopes. Moreover, we show that the relative contributions 
to the r-process abundances in both the solar-system and metal-poor stars from the weak, main, and 
fission-recycling environments required by this proposal are consistent with estimates of the relative 
Galactic event rates of core-collapse supernovae for the weak and main r-process and neutron star 
mergers for the fission-recycling r-process. 

Subject headings: nuclear reactions, nucleosynthesis, abundances - supernovae: general - stars: abun¬ 
dances 

by the relative /3-decay rates along this r-process path., 
i.e., longer /3-decay lifetimes result in higher abundances. 

In spite of this simplicity, however, the unambiguous 
identification of the site for the r-process nucleosynthe¬ 
sis has remained elusive. Parametrically, one can di¬ 
vide current models for the r-process into three scenar¬ 
ios roughly characterized by the number of neutron cap¬ 
tures per seed nucleus (n/s). This parameter, in turn 
is the consequence of a variety of conditions such as 
time-scale, baryon density, average charge per baryon. 
Ye = {ZjA), and entropy (or baryon to photon ratio) cor- 
responding to different astrophysical environments (e.g . 

IMever fc BrownlllQQTtlOtsuki. Kaiino fc Mathewsll2003ll . 

An environment in which there are few neutron cap¬ 
tures per seed (n/s ^ 50) produces what has been identi¬ 
fied a s the weak r-process (|Wasserburg. Busso fc Gallinol 
Il996f) . It can only produce the lightest r-process nuclei 
up to A ^ 125. Such an environment may occur, for 
example, in the neu trino-driven wind of core-collapse su - 
pernovae (CCSNe) (|Wooslev. et al.l[T9^IWanaio ll2013li . 
part of the o utflow from the remnant of compact bi¬ 
nary merg ers (iRosswog et al.l 12014 iPerego et al.l 12014 
lJust et akl 1201511 . the delayed magnetohydro dynamic 
(MHD) jet from CCSNe (|Nishimura et al.ll2015r i. 

An enyironment with enough neutron captures per 
seed (n/s ~ 100) to produce the two r-process abun¬ 
dance peaks at A = 130 and 195, corresponds to what 
has been dubbed the main r-process and could cor¬ 
respond, for example, to the ejection of neutron-rich 
material yia magnetic turbulence in magnetohydrody- 
namically driyen jets (MHDJ) from core collapse super- 


1. INTRODUCTION 

It has been known for more than half a century 
(jBurbidge et al.1 [l957[) that about half of the elements 
heayier than iron are produced yia rapid neutron capture 
(the r-process). Indeed, the basic physical conditions for 
the r-process are well constrained (jBurbidge et al.l[T957[) 
by simple nuclear physics. The obseryed abundance dis¬ 
tribution requires a sequence of near equilibrium rapid 
neutron captures and photoneutron emission reactions 
far on the neutron-rich side of stability. This equilib¬ 
rium is established with a maximum abundance strongly 
peaked on one or two isotopes far from stability. The rel- 
atiye abundance of r-process elements is then determined 
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novae (iNishimura et alJl2006l: iFuiimoto et al.ll2007ll2008l: 

Onoe taLl [201^ Wintel er etalJ 120121 : iNakamura et all 
2014 INishimura et alJ 2015). or in neutron star me rg- 


ers tNSMs') (jWanaio. et al.H20l4 iGorielv et al.l[2015ti . 

In this paper, we are particularly interested in a third 
environment that we dub the fission-recycling r-process. 
In this environment the number of neutron captures per 
seed nucleus can be very large {n/s ~ 1000). The r- 
process path then proceeds along the neutron drip line all 
the way to the region of fissile nuclei {A Ri 300) where the 
r-process is terminated by beta- or neutron-induced fis¬ 
sion. Fission recycling can then occur whereby the fission 
fragments continue to experience neutron captures until 
beta- or neutron-induced fission again terminates the r- 
process path. After a few cycles the abundances can 
become dominated by the fission fragment distributions 
and not as much by the beta-decay flow near the closed 
shells. Hence, a very different mass distribution can en¬ 
sue. Such environments are often associated with the 
dynamical ejecta from NSMs in which the tidal ejection 
of neutron-rich material during the merger can lead to fis- 


Gorielv. Bauswein h Jankall20: 

Lit iKorobkin et aLll2012t 

Piran. Nakar & Rosswoell20I3l: 

Rosswoff. Piran & Nakarl 


20131 iGorielv et al.ll201 


Observations ( Sneden. Cowan fc Gallinol l2008f) show¬ 
ing the appearance of heavy-element r-process abun¬ 
dances early in the history of the Galaxy seem to fa¬ 
vor the short progenitor lifetime of CCSNe over NSMs 
as the r-process site. However, identifying the r- 
process site in models of CCSNe has been difficult 
(lArnould. Gorielv fc Takahashil l2007t iThielemann et aH 


The three types of environments, neutrino-driven 
winds (NDWs), magnetohydrodynamic jets (MHDJs), 
and neutron star mergers (NSMs), have all been stud¬ 
ied extensively in the recent literature. The robustness 
of the results varies for different environments, but un¬ 
certainties in both astrophysical conditions and nuclear 
input are well recognized in all cases. For exam ple, 
the previously popular model (jWooslev. et al.l[T994[ l of 


r-process nucleosynthesis in the ND W above the newly 
forming neutron star h as been shown (|Fischer et al.l[2Mcil : 
iHiidepohl et al.l 1201011 to be inadequate as a main r- 
process site when modern neutrino transport methods 
have been employed. The required neutron captures per 
seed do not occur in the neutrino energized wind. Nev- 
erthele ss, it is quite li kely that the weak r-process does 
occur (|Wanaio 1 1201311 in the NDW producing neutron 
rich nuclei up to about A ~ 125. 

Regarding the weak r-process, however, one should 
note that while most s-process models now produce r- 
process residuals for A^120 which look remarkably sim¬ 
ilar, this is not the case for lighter nuclei. Hence, there 
is some uncertainty in determining the solar-system r- 
process abundances via subtraction of the s-process con¬ 
tribution from the total abundances. Indeed, a final 
consensus has not yet been reached on the predicted 
s-process abundances of light elements. Hence, what 
is usually taken to be the weak r-process [or unknown 
lighter element primary process (LEPP)] may actually 
corres pond to a much differe nt environment. For exam¬ 
ple in iTrippella et al.l l)2014ll it was demonstrated that 
enhanced light-element abundances could arise via non- 


parametric MHD-driven mixing mechanisms. This en¬ 
hanced light-element s-process component could obfus¬ 
cate the need for a weak r-process. Nevertheless, with 
this caveat in mind we adopt the NDW in supernovae as 
representative of the weak r-process. We note, however, 
that our arguments below regarding the relative contri¬ 
bution of the weak r-process may in fact refer to the 
relative contribution of an s-process driven LEPP. 

Indeed, the difficulties in reproducing the r- 


ies of NSMs (e.g. 

Gorielv. Bauswein Kr, Jankal 120111 

Korobkin et al.ll20I2l 

Gorielv et al.ll20I3l: IWanaio. et al.l 

20l4IPerego et al.ll20l4). Nevertheless, one scenario for 


the r-p rocess in GGSNe remains viable. It is the MHDJ 
model (INishimura et al.l20^lFuiimoto et all2007L 120081: 

OnoetlaU [ 201^~ Winteler et al. 2012t INakamura et all 

2014 12014 INishimura et al.ll201Hl . In this model mag¬ 


netic turbulence leads to the ejection of neutron rich ma¬ 
terial into a jet. As the jet transports this neutron-rich 
material away from the star it can undergo r-process nu¬ 
cleosynthesis in a way that avoids the low neutron-to- 
seed ratios associated with neutrino interactions in the 
NDW model. Moreover, the required conditions of the 
r-process environment (timescale, neutron density, tem¬ 
perature, entropy, electron fraction, etc.) are well ac¬ 
commodated in this model. 

However, there is a persist ent problem in this 
model, or any general model (e.g . iMever fc Brownlll997t 
lOtsuki. Kaiino fc Mathe^l2003 1 in which the r-process 
elements are produced in a short time scale via the rapid 
expansion of material away from a neutron star. In such 
models, the neutron density rapidly diminishes and r- 
process path freezes out near the neutron closed shells far 
from stability. Most such models underproduce isotopic 
abundances just below and above the r-process abun¬ 
dance peaks as we describe in more detail below. 

Indeed, all r-process models are fraught with uncer¬ 
tainties in the input nuclear physics, the astrophysical en¬ 
vironment, and the galactic chemical evolution. Rather 
than to give up, however, it is highly desirable to explore 
any possible method in which the relative contributions 
of each of the primary environments (weak, main, and 
fission recycling) could be ascertained from observation. 
In this paper we propose such a possibility. 

With this in mind our goal is to analyze the general 
advantages and disadvantages of each of the characteris¬ 
tic environments. Although we have noted here specific 
astrophysical models that are likely to be associated with 
the various conditions, the readers should be aware of the 
uncertainties involved and consider these environments 
as illustrative, not definitive. Nevertheless, we speculate 
here that it may be possible to quantify the relative con¬ 
tribution of each scenario to the observed solar-system 
r-process abundance distribution and the distribution of 
r-process elements in the early Galaxy. The novel con¬ 
clusion of this paper is that one can possibly utilize the 
inherent shortcomings of the three characteristic envi¬ 
ronments to estimate the relative contributions of each 
(weak, main, and fission recycling) to the final observed 
r-process abundances. 

2. EFFECT OF NUCLEAR CLOSED SHELLS 

Figure [T] illustrates why the abundances below and 
above the r-process peaks are bypassed. It shows an 
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example of a typical calculated r-process path near the 
N = 82 neutron closed shell just before freezeout when 
the neutrons are rapidly exhausted and the abundances 
begin to beta decay back to the region of stable isotopes. 
Neutron captures and photo-neutron emission proceed in 
equilibrium for nuclei with a neutron binding energy of 
about 1-2 MeV. Above and below a closed neutron shell, 
however, this r-process path shifts abruptly toward the 
closed shell from below (or away from the closed shell 
for higher nuclear masses). This shifting of the r-process 
path toward the N = 82 neutron closed shell causes iso¬ 
topes with iV = 70 — 80 (A ^110-120) to be bypassed. 
Similarly, the A =140-147 underproduction corresponds 
to the isotopes with proton closed shell Z = 50 and 
N ^90-97 {A ~140-147). These isotopes will also be 
bypassed in the beta-decay flow as is evident on Figure 

m 


We emphasize that this is not just an arti¬ 
fact of the particular mass model employed in this 
study. Nearly all models in the current litera- 
ture with a rap i d fre e zeout (including the MHDJ 
(INishimura et all 120061: i Fuiimot o et al.l 120071 120081: 


Onojrt^ 201^ jWmteleret al.l 120121 : INakamura et all 

2014HNishimura et al.ll2015ll l show this underproduction 


if the final abundances are normalized to the abundance 
peaks. Indeed, one is hard pressed to find any model for 
the main r-process (including NDW models) in which 
this underproduction does not occur. 

One can of course contrive calculations to somewhat 
fill the dips on both sides o f the second r -abun dance 
peak. Recently, for example, iLorusso et al.l (1201511 were 
able to avoid the underproduction in a schematic high- 
entropy outflow model by incorporating a summation of 
several entropies. In such a model one can fill in the 
dips similarly to the way we propose to do this by sum¬ 
ming contributions from various physical conditions. As 
another example, calculations could fill the dips by us- 
ing the ETFSl mass m odel as displayed in Fig. 7 in 
INishimura et al.l l)2006f) l. However, these models do so 
at the cost of displacing the 2nd and 3rd peaks and/or 
underproducing (or overproducing) abundances over a 
wide mass region between the second and third peaks. 
This was also a con sistent feature i n the o riginal realis¬ 
tic NDW models of lWooslev. et al.l (|1994ll . Indeed this 
effect is apparent in almost every r-process ca l culatio n 
since the 1970s (cf. review in IMathews k, War^ (11985!) '). 

We n ote, however^ that new attempts have been pre¬ 
sented (|Kratz et al.l 120141) of r-process calculations in a 
pa rameterized NDW s cenar io based upon the models 
of iFreiburghaus et ahl (|1999!) . Making use of new nu¬ 
clear masses and beta-decay r ates from the finit e-range 
droplet model FRDM (2012) (|M611er et al.ll20l3l it was 
shown that the previous discrepancies near A = 120 are 
significantly diminished compared to the s ame calcula- 
tion based upon the previous FRDM(1992) (iMoller et al.1 
Il995f) nuclear properties. Hence, one must keep in mind 
that at least some of the apparent discrepancy may be 
due to the adopted nuclear input. 

Although it has b e en speculated for s o me time (e.g. 
Wooslev. et al.l 11994 IPfeiffer et al.l 120011 : iFarouai et al.l 


2010!) that this could be due to quenching of the strength 


of the shell closure or beta-decay rates near the closed 
neutron shell, this explanation is unlikely . Recent r- 
process calculations ()Nishimura et al.ll2012!) based upon 


new measurements (|Nishimura et al.ll201l!) of beta-decay 
half lives near the A =130 r-process path have confirmed 
the absence of shell quenching effects i n the beta flow. 
More over, recently the first ever studies (|Watanabe et al.1 
1201311 of the level structure of the waiting-point nucleus 
^^®Pd {Z =46, N =82 in Fig. 1) and ^^®Pd have been 
completed. This study clearly indicates that the shell 
closure at the neutron number N =82 is fairly robust. 
Hence, there is absolutely no evidence of the hypothe¬ 
sized quenching effects in either the beta decay rates or 
nuclear masses. One must suppose that some other res¬ 
olution of this underproduction is necessary. 

One goal of this paper is, therefore, to point out that 
a solution to the underproduction of nuclei above and 
below the r-process abundance peaks can be obtained if 
one considers that a fission recycling environment (e.g. 
NSMs) has contributed to the solar-system r-process 
abundance distribution in addition to the environments 
responsible for the weak and main r-process (like CC- 
SNe). Indeed, this novel solution not only resolves this 
dilemma but can quantify the answer to the question of 
the relative contributions to the r-process abundances 
of the weak and main r-process environments (such as 
those due to CCSNe) vs. long-duration fission-recycling 
environments (such as NSMs) . 



Fig. 1.— (Color online) Illustration of the (N,Z) path of r-process 
nucleosynthesis (blue line) for nuclei with yl ~ 90 - 150 in the 
vicinity of the N = S2 neutron closed shell and Z=50 proton closed 
shell just before freezeout of the abundances in a typical main r- 
process (MHDJ) model. Black squares show the stable isotopes. 


3. FISSION RECYCLING R-PROCESS 

For this paper we highlight the possible role of fission 
recycling to account for the underproduction problem 
above and below the r-process peaks often found in mod¬ 
els for the main r-process abundances. For our purposes 
we employ a specific NSM model although we note that 
this is illustrative of any fission-recycling r-process en¬ 
vironment. Nevertheless, the most natural current site 
for such fission recycling to occur is in the NSM mod¬ 
els. The ejected matter from NSMs is very neutron-rich 
(Ye ^ 0.1). This means that the r-process path proceeds 
along the neutron drip line all the way to the onset of 
fission recycling. As noted above, after a few cycles the 
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abundances can become dominated by the fission frag¬ 
ment distributions and not as much by the beta-decay 
flow near the closed shells. Hence, a very different mass 
distribution can ensue. 

In this regard we note that a number of recent studies 
(iGorielv. Bauswein fc ,Tankall2011l: iKorobkin et al1l2012t 

Gorielv et al.l 120131 : iWanaio. et al.l 12014 Pereeo et al l 
2014fl have indicated that the r-process in NSMs can in¬ 
volve a distribution of neutron-rich environments. Such 
models can produce a final abundance pattern that is 
similar to the solar-system r-abundances. Here, even 
though we use the term NSM model, it is intended to 
refer to the portion of the ejecta that experiences fission 
recycling, while the other ejecta is similar to a NDW or 
MHDJ like model. Hence, when we refer to the NSM 
model we really mean the ejecta that experiences fission 
recycling that fills in the bypassed abundances produced 
in trajectories that produce the main r-process. 

An important point is that models including fission re¬ 
cycling effects produce a final abundance pattern that 
is re latively insensitive to the astrophysical uncertain¬ 
ties (|Korobkin et al.l[2012[) , although the total (including 
non-recycling ejecta) abundances can be sensitive to the 
detailed model. 

Nevertheless, the distribution of nuclear fission prod¬ 
ucts can affect the abundance pattern. Hence, one 
must carefully extrapolate fission fragment distribu- 
tions (FFDs) to the vicinity of the r-process path (cf. 
iMartfnez-Pinedo et al.l[200^ lErler et al.l[2?)T^ . We ar¬ 
gue that by incorporating the expected broad distribu¬ 
tion of fission fragments based upon phenomenological 
fits to observed FFDs, the effect of the neutron closed 
shells becomes smoothed out, thereby providing a means 
to fill in the isotopes bypassed in the main r-process. 

For the present study we have made use of self con¬ 
sistent /3-decay rates, /3-delayed neutron emission prob- 
abilities, and d-de layed fission probabilities taken from 
I Chiba et al.l (|2008f) . The sponta neous fi s sion r ates and 
the a-decay rates are taken from iKoural (|2004fl . In our 
r-process calculations, /3-delayed fission is the dominant 
nuclear fission mode. Hence, for the most part other fis- 
sion modes like neu tron-induced fission can be neglected 
(|Chiba et al.ll200§l . 

To generate FFDs far from stabi lity we have made 
use of a semi-empirical model (lOhta et ah. I l2007t 
iTatsuda et al.112?)?)^ iChiba et al.112008D that well repro¬ 
duces the systematics of known fission fragment distribu¬ 
tions. This model can be naturally extrapolated to the 
required heavy neutron rich isotopes of the r-process. 
As such it is a robust alternative means to predict yields 
from fission recycling. 

A key ingredient of this model is that it can account for 
FFDs that can either be single humped, bimodal or even 
trimodal. This is achieved by a weighted superposition 
of up to three Gaussian functions: 


f{A,Ap) 


E 


,— Wi exp 
V^cr 


(^) 


( 1 ) 


where A is the mass number of each fission fragment, Ap 
is that of the parent nucleus, cr is the width of the three 
Gaussian functions, and the sum is over the possible fis¬ 


sion fragment distributions, i = L,H, M, with 


Ah = 


Ar = 


and 


Am = 


q ) 

{Ap Afloss) : 

( 2 ) 

A{Ap-Nioss) , 

( 3 ) 

{Ah+Al) 

2 

( 4 ) 


The factor Wi is a weighting given by (1—Wg) for i = L,H 
and 2u}s for i = M. The quantities lOs and a are shape 
symmetry and mass-asymmetry parameters, respectively 
as defined below. Woss is the number of prompt neu¬ 
trons. 

For the present application we include the dispersion 
in the FFDs [a = 7.0) and Nioss = 2 from measured 
experiments on actinides. The adopted fission neutron 
emission is an average value for all possible fission events. 
We have run calculations in which this number A^joss is 
varied from 2 to 8 and found that the results are nearly 
indistinguishable although a very small change is found 
below A < 100 and near the valley around A = 180. The 
atomic number and neutron number of each fission frag¬ 
ment is determined by the assumption that the proton to 
neutron number ratio is the same as that of the parent 
nucleus after correcting for prompt neutron emission, i.e. 
ZplNp = Zl[N + Ni,,,l2). ^ 

We have run calculations in which Gaussian width pa¬ 
rameter a is varied from 4 to 14 and compared with the 
result with a = 7. We found that the rare-earth peak 
changes by only -1-20%,-25% so that the abundance de¬ 
creases slightly for larger a. Although the abundances 
below A < 100 and near the valley around A = 180 
increases as a increases, these changes do not change 
the overall distribution drastically, and the conclusions 
of this article are not affected by fixing a = 7. 

The quantity a in Eqs. (2) and (3) is the average mass- 
asymmetry parameter corresponding to the valley of the 
potential energy surface of the parent nucleus near the 
scission point for nuclear fissi on. This has been calcu- 
lated in the liquid drop model (iMvers fc: Swiateci^ll999tl 
with shell energy corre ctions determined (llwamoto et al.l 
119761: iSato et ^ I1979H from the two-center shell model 
in the three-dimensional shape parameter space com¬ 
prised of a, the distance between the centers of the two- 
harmonic oscillators z, and the deformation parameter 
of the fission fragments, 5. The quantity uig is deter¬ 
mined as ujs = —0.2(14 — Va — 2.0) for 14 — 14 < 2.0 
MeV, and = 0 otherwise, where 14 and 14 denote the 
potential values at symmetric and asymmetric valleys, 
respectively, at the fragment distance z corresponding to 
scission. This approximate formula is derived to account 
for the observed rapid change between asymmetric- and 
symmetric-mass distributions around ^®®Fm, i.e. 14 is 
adjusted relative to 14 to reproduce the observed mass 
distributions of the Fm isotopes with Eq. (1). 

For illustration in the present study we have car¬ 
ried out r-process simulations in the fission recy- 
cling environm e nt fr o m the NSM outflow models o f 
Korobkin et al.l (l2012f): iPiran. Nakar fc Rosswod (1201311 : 


Rosswog. Piran fc Nakarl ( 20131) . A^ an illustration of 


the main r-process we take abundances in the ejecta from 
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the MHDJ model of iNishimura et alJ (|2012r) . For the 
weak r-process we use yields from the NDW models of 
iWanaio l (IMl . 

Our adopted NSM outflow model is derived from 
3D Newtonian smooth e d-particle hydrodynamics (SPH) 
(iKorobkin et alJ 120121: iPiran Nakar fc Rosswod l2013t 
iRosswog. Piran h Nakarl 1201311 . It" gives qualitatively 
similar results to the fission recycling r-process yields 
calculated in 3D general relativistic SPH simulations 
(|Bauswein et alJ [20131 iGorielv et al.l l20l'^ and the full 
3D g eneral relativistic simulations of rWanaio. et H] 
(IMl in the heavier mass region. The details are dif¬ 
ferent for lighter masses because of the broader range of 
neutron densities and electron fractions in the particle 
trajectories in those models. 

We emphasize that the models run here can be con¬ 
sidered as generic fission recycling models. For specific 
NSM models we utilize the trajectories from the binary 
neutron star merger of two neutron stars with M = 1.0 
Mq each. Alth ough 1.0 Mfr^ is not the typical n eu- 
tron star mass (IValentim. Rangel fc HorvathI l20lHl . it 
has been shown ( Korobkinet al.ll2012H that the result¬ 
ing abundances are nearly independent of the neutron 
star masses in the binary. The hydro dynamic simula - 
tions are based upon the SPH method ofiRosswo iOpoi, 
the equation of state (EoS) of iShen et all (|1998al lbll. an 
opacity-dependent multi-flavor neutrino leakage scheme 
(jRosswog fc Liebcnddrfcil 120031 1. and Newtonian grav¬ 
ity. We use 30 available trajectories of neutron star 
merger ejecta to calculate the nucleosynthesifl The 
ejected mass from th is binary merger is ~ 0.01 Mq 
(|K orobkin et al.ll2012h . After the end of the hydrody¬ 
namic simulation at 1 5 ms) the thermody namic 

evolution can be continued (jRosswog et al.l UoidH as a 
free adiabatic expansion. 

The reason for the use of these trajectories is that 
they are publicly available and lead to robust fission re¬ 
cycling. Moreover, a main point of this paper is the 
importance of a fraction of r-process material that in¬ 
volves fission recycling. If s ome fraction of the ejects in 
the N S M calculations a s in IGorielv. Bauswein fc Jankal 
(l2011lj: IKorobkin et al.l (1201211: IGorielv et al.l ( 201311 : 
IWanaio. et al.l (|2014ll : lPerego et aLh 201411 do not involve 
fission recycling, then this paper deals with the fraction 
of material in their models that leads to fission recycling, 
while the other trajectories would be absorbed into what 
we label as other components. 

In contrast to CCSNe, baryons participating in the 
r-process constitute a large fraction of the total mass- 
energy in the NSM ejecta. Thus, the nuclear energy 
released by the r-process nuclear reactions must be in¬ 
cluded after tgn by an entropy source term, 

dS = -eth^ {niiC^/kBT + (Xi/ksT) dYi , (5) 


where a heating efficiency p arameter eth ~ 1 is intro¬ 
duced IKorobkin et al.l (|2012[ 1 to account for neutrino en¬ 
ergy losses. 

The nucleosynthesis calculations were started at a tem¬ 
perature T = 9.0 X 10® K. At this point all nuclei are in 
nuclear statistical equilibrium, and the composition is 


® Trajectories from http://compact-merger.astro su.se/ 


completely determined from the density and charge-per- 
baryon Yg of the material ejected from the neutron stars. 
At this point the material almost entirely consists of free 
neutrons plus some heavy seed nuclei with ArsTO. 

As the temperature and density decrease, however, 
the material is evolved using an updated ve r sion o f 
the nuclear network code of iTerasawa et al.l (1200111 . 
The ne utron radiative captur e rates are as summa¬ 
rized in ITerasawa et al.l (1200111 . however, for both the 
weak and main hot r-process considered here, the abun¬ 
dance patterns mainly depend on the nuclear masses 
and beta-decay rates but not on the radiative neutron- 
capture rates. This is because the system proceeds in 
(n, 7) equilibrium until a rapid freezeout of the neutron 
abundance. On the other hand, the fission-recycling 
NSM r-process considered here depends on the radia¬ 
tive neutron-capture rates because, when the tempera¬ 
ture is low, the (n, 7) and beta-decay rates (in the so- 
called ’’cold r-process”) determine the final abundances. 
In the fission recycling model adopted here, the r-process 
path terminates in a region where beta-induced fission is 
much faster than the neutron-induced fission so that the 
r-process is always terminated by beta-induced fission. 
We note, however, that this depends upon the treat- 
ment of fission barri ers and a different treatment (e.g. 
IKorobkin et al.l 1201211 can result in a different mode of 
fission termination. 

Once the r-process path fissions, we utilize the fis - 
sion fragment distributions given in lOhta et ah. I (|2007[ 1 
and also the nucl ear masses from the KTUY model 
(iKoura et al.ll200§l . The fission barriers are extracted 
from the same KTUY model. However, since the KTUY 
model treats only the symmetric fission, we adopted here 
the two-center shell model to allow more general fis¬ 
sion fragment distributions. The KTUY model has been 
shown within the GT2 theory to reproduce recent mea- 
surement s of beta-decay half-liv es of exotic neutron-rich 
isotopes (jNishimura et al.ll2011ll . In a separate forthcom¬ 
ing paper we will summarize a detailed comparison of the 
predictions of this model with known FFDs. 

We al so note that there have been numerous stud¬ 
ies (e .g. iQtsuki. Kaiino fc MathewsI l2003t I Pfeiffer et al.l 
I 2 OOIII of the sensitivity of this type of paradigm on var¬ 
ious nuclear physics parameters. However, almost all 
MHDJ (or NDW) models (without shell quenching) show 
the abundance deficiencies on either side of the closed 
r-process abundance peaks. Moreover, the NDW and 
MHDJ supernova models often involve little or no fission 
recycling. As such, they do not depend on the details of 
fission rates and fragment distributions. 

We note, however, that our NSM calculation (as 
shown below in Fig. E]) produces a different abun- 
dance pattern than that of p r evious NSM studies 
(IGorielv. Baiiswein fc Jan^l2(llll : IKorobkin et al.l[2012l : 
IGorielv et al.ll201^ . especially in the region spanning be¬ 
tween the 2nd and 3rd r-process peaks. There are two 
reasons for this difference: 1) The fission fragment distri¬ 
butions (FFDs); and 2) the number of fissioning nuclei 
contributing to fission recycling and the freezeout of the 
r-process distribution. 

Regarding the FFDs, it has often been noted (e.g . 
Gorielv. Bauswein fc Janl^ 120111: IKorobkin et al.l 1201^ 
Gorielv et al.ll2013f) that Hie elemental abundances from 
NSM calculations depend strongly on the FFD model. 
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Admittedly this is a major uncertainty in all calculations 
of fission recycling in the r-process. As noted above, our 
FFD model is based upon the KTUY model plus a two- 
center shell model to predict both symmetric and asym¬ 
metric FFDs with up to three components. As such, 
fissile nuclei in our approach can span a wide mass range 
(A=100-180) of fission fragments. This is illustrated in 
the upper panel of Fig. [2] that shows the final abundance 
distribution compared with the FFDs of three illustrative 
nuclei. 

On the ot her hand, the models of 

iKorobkin et al.1 (j2012l l are mostly based upon 

a simple two fragment distribut ion _ as in 

iPanov. Freiburghaus fc F.-K. Thielemaniil (l2001fl lor 

altern atively the prescription of iKodama fc Takahashil 
(jl975ti '). The assumption of only two fission daughter 
nuclei tends to place a large yield near the second 
r-process peak leading to a distribution that looks 
rather more like the sola r r-process a bunda nces. 
In contrast, the FFDs of iGorielv et al.l (|2013f) are 
based upon a rather sophisticated SPY revision 
( Eanebianco_et_aD [2Q1^ of the Wilkins on fission model 
( Wilkins. Steinberg fc ChasemanI Il97b[) . The main 
ingredient of this model is that the individual potential 
of each fission fragment is obtained as a function of its 
axial deformation from tabulated values. Then a Fermi 
gas state density is used to determine the main fission 
distribution. This leads to a multiple hump FFD similar 
to that considered here, but even with up to four humps. 
Although this arguably represents a more fundamental 
approach than that employed in the present work, we 
prefer the phenomenological FFD approach here as an 
alternative means to estimate fission yields far from 
stability. 

An even more important difference between the present 
work and that of previous studies is the termination of 
the r-process path and the number of fissioning nuclei 
that contribute to fission recycling and the freezeout of 
the r-process abundances. The r-process path in our 
NSM calculations proceeds rather below the fissile region 
until nuclei with A 320, whereas the r-process path in 
(jGorielv et al.l 1201311 terminates at A ^ 278 ( or for a 
maximum (Z) for ( Korobkin et al.l[2012ll l. Moreover, we 
find that only ~ 10% of the final yield comes from the 
termination of the r-process path at N = 212 and Z = 
111, while almost 90% of the A = 160 bump shown in 
Fig. [2] comes from the fission of more than 200 different 
parent nuclei mostly via beta-delayed fi s sion. This is in 
contrast to the yields of IGorielv et al.l l)2013ll that are 
almost entirely due to a few A Ri 278 fissioning nuclei 
with a characteristic four hump FFD. This is the reason 
why they obtain a solar-like r-process like distribution. 

To illustrate this point, in the lower panel of Fig. [2] 
we compare the yields of our model with a calculation in 
which we assume that the r-process path is terminated by 
symmetric fission of nuclei with A = 285. Clearly, in this 
case a solar-like distribution is obtained similar to that of 
Refs. (IGorielv. Bauswein fc ,lankall201lt IKorobkin et al.l 
120121: IGorielv et al.11201^ This highlights the impor¬ 
tance of detailed fission probabilities along the r-process 
path. 

Finally, we note that the apparent suppression of the 
the 3rd r-process peak in our final abundances relative to 
that of other works is caused by the large increase in the 
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Fig. 2. — (Color online) Illustration of the impact of fission yields 
and fission recycling on the final r-process abundances. Upper 
panel shows the relative contributions for 3 representative nuclei 
compared with the final abundance distribution. The lower panel 
shows the same final r-process yields compared with the distribu¬ 
tion that would result if fission recycling were only to occur from 
parent nuclei at the termination of the r-process path at A = 285. 

rare earth elements resulting from the FFDs of repeated 
fission recycling. 

4. RELATIVE H-PROCESS CONTRIBUTIONS 

Figure [3] shows the main result of this paper. The 
red line on Figure |3] shows the result of our fission re¬ 
cycling nucleosynthesis simulation summed over all tra¬ 
jectories of material ejected from the binary NSM model 
adopted here. This is compared with the abundances in 
the ejecta from t he main r-proces s (blue line) from the 
MHDJ model of iNishimura et al.l (1201211 , and also the 
NDW weak r-process abundances (green line) produced 
in the NDW fr om the 1.8 Mq supernova core model of 
IWanaio I (|2013ll . 

The key point of this figure is the important role that 
each process plays in producing the total abundance pat¬ 
tern of solar-sys tem r-process abundances [black dots 
with error bars ([Gorielvl 1199911 ] . The total abundance 
curve from all processes is shown as the black line on 
Figure [31 The weighting factor fFission was determined 
from a normalization to isotopes near A=145-155 for the 
fission recycling (NSM) model. The factor fweak was de¬ 
termined from a fit to light isotopes near A=100 for the 
NDW model. The MHDJ yields were normalized to the 
second r-process peak. The best fit (black) line in Figure 
|3]is for fFission = 0.16 and fweak = 4.3, or roughly 79% 
weak, vs. 18% main, and ^ 3% fission-recycling contri¬ 
butions with some uncertainty in the different models as 
noted above. Nevertheless, these estimated relative con¬ 
tributions are at least consistent with roughly estimated 
Galactic yields as described below. 

Of particular relevance to the present study is that 
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Fig. 3.— (Color online) Average final abundance patterns for 
the fission recycling environment of NSM (red line), the main r- 
process abundances from the MHDJ model (blue line) and weak 
r-process (green line) fro m the NDW. These are compared with the 
observed IIGorielvlll999ll r-process abundances in the solar system 
(black dots). The thin black line shows the sum of all contributions. 


the underproduction of nuclides above and below the 
A = 130 r-process peak shown by the blue line is nearly 
accounted for by the fission recycling (NSM) and weak 
r-process (NDW) models. The hnal NSM r-process iso¬ 
topic abundances from our adopted model for fission 
yields exhibit a very flat pattern due to several episodes 
of fission cycling. Thus, we find that fission recycling 
has the potential to resolve most of the underproduc¬ 
tion problems for the elements just below and above the 
abundance peaks in models of the main r-process. The 
remaining underproduction below the A = 130 peak is 
most likely due to the weak r-process as illustrated on 
Figure [3l 

The main point of this paper is that one can deduce the 
relative contributions of each r-process model based upon 
their relative shortcomings. However, it is important to 
ask whether the inferred fractions, of ^79% NDW, ~18% 
MHDJ, and ~3% NSM are plausible. 

Although there are many uncertainties in the as- 
trophysical and ga l actic chemical evol ution parameters 
(|Argast et al.ll200(l : iKomiva et al.ll201^ , it is worthwhile 
to estimate weight parameters fFission and fweak from 
observed Galactic event rates and expected yields. In 
particular we write 


and 


Jfz. 


RNSklMr^NSM 

eMffnjR.CCSNMi.,MHDJ 


fweak 


RcCSNMr,Weak 

eM/fDjRcCSNMr^MHDJ 


( 6 ) 

(7) 


where Mr,NSM, Mr,MHDj, and Mr,weak are the ejected 
masses of r-elements from the NSM, MHDJ, and NDW r- 
process models, respectively, while Rccsn and Rnsm are 
the corresponding relative Galactic event rates of GGSNe 
and NSMs. 

The ejected ma ss of r-process elements in the models 
of lWanai^ 1)20131) is « 2 x 10“® Mq and nearly indepen¬ 
dent of assumed core mass. The quantity cmhdj is the 
fraction of GGSNe that result in magneto-rotat i onally 
driven jets. This was estimated in iWinteler et ahl (|2012l) 


to be ^ 1 % of the core-collapse supernova rate based 
upon the models of iWooslev &: HegCT (I2006D . However 
this is probably uncertain by at least a factor of two. In¬ 
deed, the fraction could be larger as most massive stars 
are fast rotators and the conservation of magnetic flux 
should often lead to high magnetic fields in the newly 
formed proto-neutron star. Hence, this fraction could 
easily range from ^ 1 to ~ 5% which incorporates the 
~ 1 % fraction of observed magnetars compared to nor¬ 
mal neutron stars. [We treat this as a lower limit because 
some fraction of observed normal neutron stars may have 
had a larger magnetic field in the past.] The mass of 
synthesized r-proc ess elements from MH DJs is estimated 
to be 6 xlO“^M 0 ([Winteler et al.l[ 2012 D while that of a 
typical binary NSM is expected to be 2 ± 1 x 10~^Mq 
()Korobkin et al.1 12012D . If the Galactic neutron star 
merger rate is SOtyo^ Myr“^ (iKaloeera et al.l[20M) . and 
the Galactic supe rnova rate is, 1.9 ± 1.1 x 10^ Myr“^ 
(jPiehl et al.l l 200 fil) . then one should expect /Fission ~ 
0.6 ± 0.4 and fweak ~ 8 ± 6 corresponding to relative 
contributions of ^ 80% weak, ^ 10 % main and ~ 10 
% fission recycling. Thus, although there are large un¬ 
certainties, these fractions are plausibly consistent with 
our fit parameters. This suggests that such a fit may be 
a way of constraining the relative contribution of NSMs 
and GGSNe to solar-system material. 

We note, however, that other NSM calculations pre¬ 
dict about 10 ~ ^ to 10 ~^ Mfn of r-process material to 
be ejected (e.g. iHotokezaka. et al1l2013t iBauswein et al.1 
nnn). Adopting a value of 10 ^ Mq could lead to 
/Fission ~ 0 . 02 , i.e. about an order of magnitude be¬ 
low that suggested in our fit to Figure [3] 


Of course, this needs to be better quantified in 


more detailed chemical evolution (IGescutti & Ghianoini 

2014 

Gescutti et al.1 

2015t iTsuiimoto & Shigevama 

2m4ayb 

: Komiva et al. 

12014 llshimaru et al.1 120151 

Wehmever et al. 1201 51) 

and chemodynamical stud- 

ies (Shell etalJ 120151: Ivan de Voort et alJ I2015D 


along with better r-process hydrodynamic models 
(iWinteler et al.ll2012l: iPerego et al.ll20l4 iRosswog et al. 

20141 : iWanaio. et al.ll20l4lGorielv et al.ll2015l : lJust et al. 

2015t iMshimura et al.ll2015h . Nevertheless, based upon 

the models adopted here, the inferred division of 
r-process contributions remains at least plausible. 

5. UNIVERSALITY OF R-PROCESS ELEMENTAL 
ABUNDANCES 

In the above we have not discussed a very impor¬ 
tant clue to the origin of r-process abundances. It is 
by now well established (|Sneden. Gowan fc Gallin^l2008D 
that the elemental abundances in many metal-poor stars 
show a pattern that is very similar to that of the solar- 
system r-process distribution, particularly in the range 
of 55 < Z < 70. This however can pose a diffic ulty 
((Mathews. Bazan fc Gowanlll992l : iArgast et al.ll2000t) for 

NSM models (either in the present work or in other stud¬ 
ies). That is because metal-poor stars are thought to 
have arrived very early in the history of the Galaxy, 
whereas NSMs require a relatively long gravitational ra¬ 
diation orbit decay time prior to merger (^ 0.1 Gyr). 
Whatever the situation, it is of value to examine the im¬ 
pact of the possible late arrival of fission recycling ma¬ 
terial on the r-process elemental abundance distribution 
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Fig. 4. — (Color online) Average final elemental abundances for 
the total sum from Fig. |2] (solid line) and the contribution with¬ 
out NSMs (dashed line). These are compared with the observed 
elemental r-process abundances in two well-studie d metal-poor r- 
proces s enhanced stars, HD1601617 (filled circles: IRoederer et al.l 
li2012l ii and CS22892-052 (open squares; ISneden eTlirT * I200H1 1. 
The curves are arbitrarily normalized at europium (Z=63). 

in metal-poor stars. 

Figure S] shows the elemental abundance distribution 
calculated in two scenarios, i.e. with and without the fis¬ 
sion recycling yields of NSMs. These are compared with 
the observed elemental r-process abundances in two well- 
studied metal-poor r-p rocess enhanced sta rs, HD1601617 
(IRoederer et al.l 120121) and CS22892-052 (|Sneden et ah. I 
120031) . Here, we note that there is little distinction be¬ 
tween the two curves (although the fit is slightly bet¬ 
ter when the fission recycling yields are included). The 
reason for this insensitivity is that the fission recycling 
environment only contributes about 3% to the total r- 
process abundance. Although this yield is important to 
fill in the isotopic abundances above and below the r- 
process peaks, and also to make the rare-earth bump near 
A=160, there is little apparent difference in the elemen¬ 
tal abundances with or without neutron star mergers. 
Among other things, this is because the region below the 
peak [Z ^ 50) is poorly sampled, and moreover, sum¬ 
ming over isotopes to produce elemental averages some¬ 
what washes out the underproduction above and below 
the r-process mass peaks. Hence, the elemental r-process 
abundances in metal poor stars do not clearly require 
that fission recycling occurred early in the Galaxy. 

We do note, however that the dispersion in the stars 
themselves for the lightest elements (30 < Z < 50) is 
suggestive that not all CCSNe contribute both a weak 
and main r-process. This is consistent with the expec¬ 
tation that the NDW could occur in all CCSNe while 
the main r-process from the MHDJ will only occur in a 
limited fraction of CCSNe, i.e. those with rotation and 
strong magnetic fields. 

6. DISCUSSION 

The fits to the abundance distribution (e.g. Figure [3]) 
are as good as or better than most models in the litera¬ 
ture. Nevertheless, it is worthwhile, to address some of 
the detailed deficiencies in both Figures [3] and O For ex¬ 
ample, although the r-process peaks at A=130 and 195 
along with the rare-earth peak region A = 145-180 in 
Figure [3] are remarkably well reproduced, there are some 


differences just above the main r-process peaks in the 
regions of A=140-145 and 200-205. We note, however, 
that these isotopes have the largest uncertainties in the 
r-process abundances themselves as is visible on FigurejS] 
Hence, these discrepancies may simply reflect the abun¬ 
dance uncertainties, although the possibility remains of 
a shortcoming in the models for these isotopes. 

Similarly, in Figured there is an underproduction of el¬ 
ements at Z=58 and 60. The abundance of Ce (Z=58) in 
Figure 4 is well determined observationall y for CS22892- 
052 a s follows: loge(Ce) = —0.50 ± 0.07 (Sneden et al~1 
[Mil and = -0.38 ± 0.08 (IHonda et al.ll2004l) . This 
corresponds to the deficient isotopes with A=140 and 
142 in Figure [3l However, the odd elements with Z= 
57 (La) and Z=59 (Pr) are reproduced. This suggests 
that the odd-even effect in the region of lanthanide el¬ 
ements may be underestimated in the mass model em¬ 
ployed here. Nevertheless, the main point of this paper 
is not to give a precise reproduction of r-process elemen¬ 
tal abundances but rather to demonstrate the possibility 
that fission recycling supplements the underproduced el¬ 
ements. Clearly, a better understanding of the nuclear 
uncertainties within this context is still needed. 

We also note that there is a possible deficiency of Pb 
(Z=82) in Figure m This, however, may relate to obser¬ 
vational uncertainties. Ther e are two measured Pb abun¬ 
dances for CS22892-052 in iSneden et ah. I ([20031 ). One 
was a ground-based measurement, while the other was 
obtained with HST. However, bo th of these values sho uld 
be considered upper limits. In ISneden et ah. I (I2003D it 
was noted that the suggested detections of the two Pb 
I lines in the ground-based spectra should be nearly 10 
times weaker than the A = 2833.05 line, that could not 
be detected in the HST spectrum. Hence, the derived 
Pb abundance upper limit from the A = 2833 line is 
probably more reliable than the abundances determined 
from the questionable detections of the other two Pb 
I lines. Thus, one should abandon the Pb abundance 
of loge(Pb) = 0.05 from the ground-based observation 
in favor of loge(Pb) < —0.2 from the HST observa- 
tion. We also note th at the more recent observation of 
IRoederer et al.l (|2009f) also obtains loge(Pb) < —0.15. 
These upper limits are consistent with our calculation. 

Another issue worthy of discussion is that of Th (Z=90) 
and U (Z=92) production in Figure [4l Th has been ob¬ 
served in a number of metal-poor stars and U in a few. 
This indicates that the r-process mechanism at work in 
the early Galaxy could produce the actinide elements 
and beyond. Although one tends to think that the pro¬ 
duction of actinide elements requires a fission-recycling 
r-process, in fact the production of Th and U is possi¬ 
ble even in models that do not lead to fission recycling. 
For exa mple, the MHDJ mode ls with strong magnetic 
fields in iNishimura et al.l (|2015[ ) could produce Th and 
U in as much as their solar abundances. On the other 
hand, MHDJ models with weak magnetic fields tend to 
produce actinides below solar abundances. Hence, the 
observation of Th and U in metal-poor stars constrains 
the early astrophysical environment, but does not neces¬ 
sarily require that a fission-recycling r-process (such as 
the NSM model) contributed to metal-poor stars. 

7. CONCLUSIONS 
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In summary, we have considered the relative contribu¬ 
tions of three generic r-process environments to the solar- 
system r-process abundances and the abundances in r- 
process enhanced metal-poor stars. These environments 
are discussed in the context of neutron star mergers, neu¬ 
trino driven winds and magnetohydrodynamically driven 
jets, although these environments should be considered 
as illustrative and not definitive of the specific r-process 
environments. Nevertheless, based upon our adopted fis¬ 
sion fragment distributions we find that the relative con¬ 
tributions from each environment has the possibility of 
explaining a unique feature of the r-process abundances. 
Moreover, the deduced relative contributions are plausi¬ 
bly consistent with galactic chemical evolution consider¬ 
ations. 

Clearly, more work along this line is required to ex¬ 
plain details. Nevertheless, we suggest that the possi¬ 


bility that all three general environments occur in de¬ 
tectable amounts in the r-process distribution should be 
taken seriously in future investigations of the origin of 
r-process nuclides. 
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